Benchmark of Schemes for Multiscale Molecular Dynamics Simulations.

In multiscale molecular dynamics simulations the accuracy of detailed models is combined with the efficiency of a reduced representation. For several applications - namely those of sampling enhancement - it is desirable to combine fine-grained (FG) and coarse-grained (CG) approaches into a single hybrid approach with an adjustable mixing parameter. We present a benchmark of three algorithms that use a mixing of the two representation layers using a Lagrangian formalism. The three algorithms use three different approaches for keeping the particles at the FG level of representation together: 1) addition of forces, 2) mass scaling, and 3) temperature scaling. The benchmark is applied to liquid hexadecane and includes an evaluation of the average configurational entropy of the FG and CG subsystems. The temperature-scaling scheme achieved a 3-fold sampling speedup with little deviation of FG properties. The addition-of-forces scheme kept FG properties the best but provided little sampling speedup. The mass-scaling scheme yielded a 5-fold speedup but deviated the most from FG properties.

Stochastic Dynamics with Correct Sampling for Constrained Systems.

In this paper we discuss thermostatting using stochastic methods for molecular simulations where constraints are present. For so-called impulsive thermostats, like the Andersen thermostat, the equilibrium temperature will differ significantly from the imposed temperature when a limited number of particles are picked and constraints are applied. We analyze this problem and give two rigorous solutions for it. A correct general treatment of impulsive stochastic thermostatting, including pairwise dissipative particle dynamics and stochastic forcing in the presence of constraints, is given and it is shown that the constrained canonical distribution is sampled rigorously. We discuss implementation issues such as second order Trotter expansions. The method is shown to rigorously maintain the correct temperature for the case of extended simple point charge (SPC/E) water simulations.

Efficient Algorithms for Langevin and DPD Dynamics.

In this article, we present several algorithms for stochastic dynamics, including Langevin dynamics and different variants of Dissipative Particle Dynamics (DPD), applicable to systems with or without constraints. The algorithms are based on the impulsive application of friction and noise, thus avoiding the computational complexity of algorithms that apply continuous friction and noise. Simulation results on thermostat strength and diffusion properties for ideal gas, coarse-grained (MARTINI) water, and constrained atomic (SPC/E) water systems are discussed. We show that the measured thermal relaxation rates agree well with theoretical predictions. The influence of various parameters on the diffusion coefficient is discussed.

Surface UV aging of elastomers investigated with microscopic resolution by single-sided NMR.

Depth profiles taken from the surface of UV irradiated natural rubber sheets have been measured with microscopic resolution using a Profile NMR-MOUSE. An NMR observable related to the sum of the spin echoes in the Carr-Purcell-Meiboom-Gill pulse sequence was used to characterize the cross-link density changes produced by the action of UV radiation in each sheet. The aging process was investigated as function of irradiation time and penetration depth. An exponential attenuation law with a space dependent absorption coefficient describes the change in the NMR observable with penetration depth. An Avrami model is used to describe the dependence of the absorption coefficient on the aging time. The method can be applied to investigate the effect of various aging agents on the surfaces of elastomers.